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The model is capable of approximating most of the Known types of 


rate-independent symmetrical hysteresis loops encountered in the 
practice. The model allows building smooth, piecewise-linear, 
hybrid, minor, mirror-reflected, inverse, reverse, double and triple 
hysteresis loops. The error of approximation of a hysteresis loop by 
the model does not exceed 1% as a rule 
































With the original model, a family of hysteresis loops is described by 
the following parametric equations 


x(a) = acos” a +b,sin" a, 
y(a)=b, sina, 


(1) 


where ais a real parameter (a=0...27); ais x coordinate of the split 
point; 6,, b, are the saturation point coordinates; m is a positive odd 
integer (m=1, 3, 5, ...) defining the curvature of the hysteresis loop; 
nis a positive integer defining the type of the hysteresis loop and its 
curvature. With n=1, the Leaf loop type is formed; with n=2, 4, 6, ... 


— the Crescent (Boomerang), and with n=3, 5, 7, ... — the Classical. 





Hysteresis loops of Leaf 
(n=1), Crescent (Boomerang, 
n=2), and Classical (n=3) 
types. The area of all the 
three loops Is the same 











Hysteresis loop (1) can always be represented as a sum of two 
parametric curves 


x(a) = X, (a) + X5 (a), 
y(a) = y,(ar) + yp (ar) 


where x,(@)=6,sin"@, y,(a)=b,sina@ is the unsplit loop; 
X,(@)=acos'a, Y,(a)=0 Is the splitting curve. 














The generating function x(q) can easily be expanded in the Fourier 
series 


x(a) =¥(A, cos(ka) +B, sin(ka)), 


ka (3) 
y(a) =b, sina, 


where the Fourier coefficients A,, B, for odd n are determined by the 
ger formulas 

A, = amt 4 oe? 

B, = (1) "2 oe O,8 2 

where Cy = Ml/[kI(/- ¥) is a binomial coefficient. 


(4) 





Having the Fourier coefficients A,, B,, the generating function x(a) 
can also be represented as 


x(a) = Sy Am, cos(ka — 9, }, (5) 


k=1 


where amplitudes Am, and phases @, of the harmonics are 
determined by the following formulas 
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The required loop tilting by angle @at the split point a is performed 
using the following transformations 


x(a) = x(a) + sin6(b, sind + b, cos o\sina—sin’ a) 7) 
(a) = y(@)+ sin@(b, cos @ — b, sin o\sina—sin’ a) 


Area S of the loop (7) is calculated by the formula 
n-1 


2 
S= aoe sin 6(b, cos@—b,sin@ | im+ NIT aie 
k=0 


a m+n-—2k 





Phase shitts 


Phase shifts Aa@,, Aa, Aa, permit us to tilt the hysteresis loop 
smoothly by the required angle @at the split point a as well as to 
smoothly change the curvature of the loop 


x(a) = 4cos”(a+ Aa@,) + b, sin"(a@+ Aa,), 
y(a) = b, sin(a+ Aa), 


where a, b, are corrected parameters of a, b,, respectively. 


The corrected parameters 2, b, can be found by the formulas 
acos”(Aa, — Aa,)—b,sin’ (Aa, — Aa,) 
sin” (Aa, — Aa, )sin"(A@, — Aa, )+cos”(Aa, — Aa, )cos"(Aa, — Aa,)’ 
5 - asin” (Aa, — Aa,)+b,cos” (Aa, — Aa, ) 
* sin” (Aa, — Aa, )sin’(Aa@, — Aa, )+cos”(Aa, — Aa, )cos"(Aa, — Aa, ) 


-o 
(10) 














1 x(a) 


(Cc) 


(a) Tilting with phase shift A@,, continuous change in the curvature by 
phase shift (b) Aa, (c) Aa, 
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>= 1) by OF 6 
—C,,2 cos(kAa@,)+(-1) 2 a C_? sin(kAa, ) 
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— Ea Boe 
C2 sin(kAa,)+(—1) 2 - C2 cos(kAa,) 
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A m-1 a n-1 


_ C2 cos(Aa, — Aa, )+ an C,? sin(Aa, — Aa, ) mb, 


>= 


_ Am,(cos Aq, — tang, sinAa, ) 


tan’ @, +1 / 


= (A.cos Aa, — B,sinAa, )ab, 


70D 


Note that the loop area depends on the amplitude and phase of the 
1st harmonic only; the rest of harmonics do not affect the loop area. 





Hysteresis loops can be tilted by skewing of the coordinate system by 
angle @along the x axis and by angle xalong the y axis 


X(ar) = x(a) + tano(b, tan« + b, \sin@ — sin” a) 


13 
y(a)= y(a@)+ b, tanx(sina — sin’ a) (13) 


The area of the skewed loop is calculated by the formula 


m_1 


2° 
a= 74 Co? b. t 1- 1)} | ————_ 
me aan (m+ 1) eee 



























































Tilting the Classical loop Changing the curvature of the Classical loop 

by skewing the by skewing the coordinate system by angle «x 
coordinate system by along the y axis. Tilting the loop at the split 

14 angle @along the x axis point by angle (a) G@15°, (b) G-15° 











Replacing the sine and the cosine in the generating functions x(q), 
y(@) of model (1) with unit-amplitude trapezoidal pulses trp,, trp,, 
respectively, one can produce piecewise-linear hysteresis loops built 
on trapezoidal pulses 


x(a) = atrp” a+(b, — a)trp” a, 
(15) 
y(@) = by, trp, Q, 


where a=0.../; T is the pulse period. 








The trapezoidal pulses tro are defined as follows 


i 


r 
trp, @ = pa + ) 


trp. a= YI 5*al2 = rit 1) rect, (@,i)+ (—1)' rect, (a, ), 


where d, D are the upper and the lower bases of the trapezoidal 


pulses, respectively (D=3d, T=d+D=40); rect,, rect. are rectangular 
pulses. 














The rectangular pulses rect,, rect, are determined by a step function 


H(@) (Heaviside function) 


rect,(a@,i)= Ha-Fr " — *)- Ha 


rect, (a, i) = Ha-Fi-278)- Ha-Fr oe) 
2 4 2 4 


























Piecewise-linear hysteresis 
loop Leaf (Play without 
Whiskers, n=1), hybrid 

Crescent (hybrid Boomerang, 
n=2), and hybrid Classical 
(n=3) built on trapezoidal 
pulses. The area of all the 

three loops Is the same 





hase shifts A@ 


Taking into account the phase shifts A@,, Aa, Aa, equations (15) 
are written as 


x(a) = atrp" (a + Aa,)+ b, trp’ (a+ Aa, ), 18) 
y(a) = b, trp,(@ + Aag). 


The corrected parameters 2, b, are determined by the following 
formulas 


- atro’(a@, + Aa, — Aa@,)—b, trp? (Aa, — Aa, ) 
~ trp” (Aa, — Aa, )trp” (a, + Aa, — Aa, )— trp” (a, + Aa, — Aq, )trp"(Aa@, — Aa, )’ 
- b,trp” (Aa, — Aa,)— atrp” (a, + Aa, - Aa,) (19) 
“trp” (Aa, — Aa, )trp2(a, + Aa, — Aa,)— tro” (a, + Aa, — Aa, )trp’ (Aa, — Aa) 

















To obtain piecewise-linear loops with gain/attenuation % one should 
add to (18) an extra term (curve) responsible for gain/attenuation 
(m=n=1) 


x(a) = x(a), 
y(a) = y(ar) + tany|x(@)— b, trp, (a+ Aa, )| 


(20) 
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Piecewise-linear hysteresis loops with gain/attenuation xv. (a) Play-Relay-Play, 
(bo) Play-Play, (c) Play-Relay built on trapezoidal pulses using the phase shitts 

















Whiskerless Play and Non-ideal Relay loops with gain/attenuation are 
built by the following equations 


(a) = atan Gtrp, a+ (b, — b, tan y)trp. OL 
- tan 6 —tany 
y(a) = tan B(x(a)— atrp, a) 


J 


where the loop parameters a, b,, b,, Pare interrelated as 
tanf=b,/(b,-a). 














The areas of the simplest piecewise-linear loops (21) having a 
parallelogram shape are determined by the formula 


$= 44 


(a— b,)tany+ by, B- asin #siny 
tan B-tany sin(B—- 7) 





Whiskerless hybrid Classical loops with the desired slope f=1/2- 6 at 
the split point are built according to the following transformation 


X(a) = x(a)+ by sin a(tro. @—trp” @) 


- (23) 
y(az) = y(@), 

Whiskerless hybrid Classical loops having the required slope fat the 
Split point, the required inclination (gain/attenuation) vot the 
rectilinear section, and the required curvature «of the curvilinear 
section are built according to the following transformation 

x(a) = x(a) +tano|(b, - a)tan«— b, tany+ b (trp, a@- tro’ a} 

V(a)= y(a)+[(b, — a)tan— b, tany|trp. a— trp’ a)+ atany(trp, atrp’ a— trp’ a) 
where k=2, 4, ... Is an additional parameter to control the loop 
curvature. 


24) 
























Hybrid Classical loops Hybrid Classical loops with specified slope 
with the specified slope b=n/2-0, gain/attenuation % and curvature « (a) 


B=1/2-0 at the split point Various gains yvfor fixed Gand x, (bd) various 
curvatures «for fixed Gand v 


26 


The area of hybrid Classical loop (23) is calculated by the formula 


S = 4ab,. (29) 


The area of hybrid Classical loop with non-zero gain/attenuation (24) 
is calculated by the formula 


k(n—1) 


(k +1)(n +k) 


+(b, ~ al K tan y —tan ‘| 
n+k 


\(b, — a)tan«— b, tany + b | 


tanétany + 1 


S=4a (26) 





Loops built on triangular pulses 


Besides trapezoidal pulses tro, formulas (9) can operate with 
triangular pulses tri, which are particular cases of the trapezoidal 


pulses (Q=0, /=D) 

x(a) = atri” (w+ Aa,) + b, tri.(@+ Aa,), 

y(a) = b, tri,(a+ Aa). 
The corrected parameters 2, b, are determined by the following 
formulas 

atri?(Aa, — Aa, )—b, tri? (Aa, — Aa,) 

tri” (Aa, — Aa, )tri?(A@, — A@,) + tri” (Aa, — Aa, )tri?(Aa@, — Aa, ) 
6 - atri” (Aa, —Aa,)+b,tri” (Aq, — Aa, ) 

“tri? (Aa, — Aq, )tri?(Aa, — Aa, )+ tri” (Aa, — Aa, )tri? (Aa, — Aa, ) 


a= 


(28) 




















The triangular pulses tri are defined as follows 


/ 


tri, @= Fy (e-F!|e- 1) rect(a, i) 


tri, @= i,(a + r} 
4 


where rect are rectangular pulses. 
The rectangular pulses rect are determined by a step function H(q) 


\oul gt ja fl lH gL j_t 
ect(ani)=H a cist H(a a r) (30) 

















Piecewise-linear loop Leaf 
(bi-linear loop, Play without 
Whiskers) at different values 
of phase shift Aa, built on 
triangular pulses 
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The general expression describing a piecewise-linear hysteresis loop 
Play with Gain is as follows 

x(a) =b, tri. a, 

a, tan B (31) 


——4_—_*__ |+) tanytri, @, 
tan 6-tany 


y(a) = (b,—b,tan /he.| a 7 


where @,=a!1/(4b,) is the value of parameter q@at split point a. 
The upper base d of the trapezoidal pulses trp, in (31) is determined 
according to the formula 
— T(b,tanB-b, ) 
- 2b, (tanB—-tany)’ 
and the lower base D according to the formula D= /-d. 


(32) 

























(a) Play with 
Gain, (b) Play 
with Gain w/o 
Whiskers 
(parallelogram 
loop), (c) Play 
with Attenuation, 
(d) Play with 
Attenuation w/o 
Whiskers, (e) 
Play (backlash), 
(ft) Play w/o 
Whiskers 









































(g) Non-ideal 
Relay with Gain, 
(hn) Non-ideal 
Relay with Gain 
w/o Whiskers, (1) 
Non-ideal Relay 
with Attenuation, 
“ (Y) Non-ideal Relay 
with Attenuation 
w/o Whiskers, (k) 
Non-ideal Relay 
(Schmitt trigger), 
(1) Non-ideal Relay 
w/o Whiskers 
_| (rectangular loop) 


The equations for a double smooth loop non-self-crossing in the 
origin of coordinates (O-shaped loop) are as follows (q=0...27) 


X(a) = {2a 5 san(e it) Nas) + b, sgn(z — a) 


= (2recta- iL {20 —~ Aa — 4 + b.} 


y(a@) = y(2a 7 5 San(e ~a)- Nas) + b, sgn(z — @) 


= (2recta— iL y(2a —~ Aa — 4 ~ b} 


where sgna=a/|aq| is the signum function; recta@=H(a@)-H(a-7) is 
a m-wide rectangular pulse. 











The loop with a self-crossing in the origin of coordinates (8-shaped 
loop) can be built according to the formulas (a=0...27) 


X(a@) = X (24 - 4 sgn(z — a)— say + b, sgn(z — a) 
= rect oo 20 — AQ, — 4 + b + (1—rect a (3 — AQ, — 2a = b.} 
y(a) = y (20 ~ 4 sgn(z — a)— say + b, sgn(z — @) 


= rect of y(2a —AQ@, — 4 + by + (1-rect a 2 —Aa@, — 2a — b} 


The loops (83), (84) have no appearance differences. By replacing x with 
T/2 (a=0... 7), double piecewise-linear loops are built by formulas (33), (34). 


(34) 





ible hysteresis loops 














2 1.5 “ 15  2X(a) 2 1.5 a 15  2X(a) 


(b) 


Double (a) smooth, (b) piecewise-linear hysteresis loop formed by linking 
36 two loops (a) Classical, (bo) Play with Gain in the saturation point b 








a=0 
aq = 0.14 (blue) 


Double non-seltf-crossing 
hysteresis loop of Propeller 
type formed by pinching a 

loop at the point of 
coordinate origin by means 
of zero splitting aanda 
phase shift Aa, 





A triple loop is assembled of three loops — one central loop and two 
outside loops linked at the saturation points b (a=0...27) 


X(@) = (recta +rect(a— z))x. (3a ~ 4 


+ (recta +rect(a —z)- i) x (2 3a - 4 —(b,, + b,,.)sgn(z — 2) 


(a) = (recta +rect(a—z))y, (3a _ 4 


+ (recta +rect(a —z)- J AE 8a - | —(b,, + b,, )sgn(z - 2 


where x,(a@), y,(@) are equations of the central loop; x,(@), y.(a@) are 
equations of the outside loops; 0,,, 6,,, are coordinates of the saturation 
point of the central loop; 6,,, 64, are coordinates of the saturation points 
of the outside loops; recta=H(a@)-H(a-m/3) is a m/3-wide rectangular pulse. 











When assembling loops (35), the condition 4= % is usually met, where 
¥%,, % are slope angles of tangents to the unsplit (a,=0) central loop and 
the unsplit (a,=0) outside loops, respectively, at the saturation point 0,. 
The slope angles 7, % of the tangents are defined by the formulas 


b 
y, = arctan—“, 
N, 1x 
(36) 


y= amin. 
Ny 2X 
In case the argument 3q@ot functions x.(@), Yo(@) in equations (35) is used 
with the plus sign, a non-self-crossing loop is obtained, and in the case of 
the minus sign — a self-crossing loop; both the loops have the same 
appearance. Triple non-self-crossing and self-crossing piecewise-linear 
loops are built by formulas (35) by replacing xz with //2 (a=0... 1). 



































































Triple (a) smooth, (b) piecewise-linear hysteresis loop formed by linking three 
loops (a) Classical, (b) Play with Gain in the saturation points b. In the linking 
points, the loop can be made both non-self-crossing and self-crossing 


Triple hysteresis loops (35) are useful for producing single smooth or 
single hybrid hysteresis loops having long whiskers. The triple loop 
consists of the central Tilted Classical loop and whiskers, which are 
formed from a pair of the outside unsplit loops of Leaf type oriented at 
the angle %. To obtain hybrid Classical loop with whiskers, the 
following equations can be used (a=0... 1) 


= (rao -Z) (20-2 
Deed) mstet] 
he 


+ 
2 
1) v.(30- | — (b,, + b)san| 2 — a 


r 
+|rect@a+rect| a - 5 


+|rect@a+rect]| ~- 





paso 
Ft =lroe rece 
(emia 






























simulation of a single Classical (a) smooth, (b) hybrid loop with long whiskers 
by means of a triple non-self-crossing hysteresis loop. Whiskers are the 
outside pair of (a) smooth, (b) piecewise-linear unsplit loops of Leaf type 


Squeeze causing a fo 


The triple self-crossing loops 
are formed by setting up a 
negative phase shift Aq, (or 
positive Aa.) that “squeezes’ 
the loop so tight that a 
foldover appears 
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